Effect of nonlinear filters on detrended fluctuation analysis 
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When investigating the dynamical properties of complex multiple-component physical and physi- 
ological systems, it is often the case that the measurable system's output does not directly represent 
the quantity we want to probe in order to understand the underlying mechanisms. Instead, the 
output signal is often a linear or nonlinear function of the quantity of interest. Here, we investigate 
how various linear and nonlinear transformations affect the correlation and scaling properties of a 
signal, using the detrended fluctuation analysis (DFA) which has been shown to accurately quan- 
tify power-law correlations in nonstationary signals. Specifically, we study the effect of three types 
of transforms: (i) linear (j/i = axi + b); (ii) nonlinear polynomial (j/i = ax^); and (iii) nonlinear 
logarithmic [yi — log{xi + A)] filters. We compare the correlation and scaling properties of signals 
before and after the transform. We find that linear filters do not change the correlation properties, 
while the effect of nonlinear polynomial and logarithmic filters strongly depends on (a) the strength 
of correlations in the original signal, (b) the power k of the polynomial filter, and (c) the offset A 
in the logarithmic filter. We further apply the DFA method to investigate the "apparent" scaling of 
three analytic functions: (i) exponential [exp(±a; -I- a)], (ii) logarithmic [log(a:: -I- a)], and (iii) power 
law [{x + a)^], which are often encountered as trends in physical and biological processes. While 
these three functions have different characteristics, we find that there is a broad range of values for 
parameter a common for all three functions, where the slope of the DFA curves is identical. We 
further note that the DFA results obtained for a class of other analytic functions can be reduced 
to these three typical cases. We systematically test the performance of the DFA method when 
estimating long-range power-law correlations in the output signals for different parameter values in 
the three types of filters and the three analytic functions we consider. 

PACS numbers: 05.40.-a 



I. INTRODUCTION 

Many physical and biological systems under multi- 
component control mechanisms exhibit scale-invariant 
features characterized by long-range power-law correla- 
tions in their output. These scaling features are often 
difficult to quantify due to the presence of erratic fluctu- 
ations, heterogeneity, and nonstationarity embedded in 
the output signals. This problem becomes even more 
difficult in certain cases: (i) when we cannot probe di- 
rectly the quantity of interest in experimental settings, 
i.e., the measurable output signal is a linear or nonlin- 
ear function of the quantity of interest; (ii) when mea- 
suring devices impose a linear or nonlinear filter on the 
system's output; (iii) when we are interested not in the 
output signal but in a specific component of it, which is 
obtained through a nonlinear transform (e.g., the mag- 
nitude or the sign of the fluctuations in the signal); (iv) 
when comparing the dynamics of different systems by ap- 
plying nonlinear transforms to their output signals; or (v) 
when pre-processing the output signal by means of linear 
or nonlinear filters before the actual analysis. Thus, to 
understand the intrinsic dynamics of a system, in such 
cases it is important to correctly analyze and interpret 
the dynamical patterns in the system's output. 

Conventional two-point correlation, power spectrum, 
and Hurst analysis methods are not suited for nonsta- 
tionary signals, the statistical properties of which change 



with time |ll |2j 13 • To address this problem, detrended 
fiuctuation analysis (DFA) method was developed to ac- 
curately quantify long-range correlations embedded in a 
nonstationary time series 0, Q- This method provides 
a single quantitative parameter — the scaling exponent 
a — to quantify the scale-invariant properties of a sig- 
nal. One advantage of the DFA method is that it al- 
lows the detection of long-range power-law correlations 
in noisy signals with embedded polynomial trends that 
can mask the true correlations in the fluctuations of a sig- 
nal. Recent comparative studies have demonstrated that 
the DFA method outperforms conventional techniques in 
accurately quantifying correlation properties over a wide 
range of scales jSj iS: iS S EyL, The DFA method has 
been widely applied to DNA^R liJ3jn|jgjllL car- 
diac dynamics ra El El El H lilllMililll H 

I25I 129 . I27I I28I 129 . ISOll. human electroencephalographic 
(EEG) fluctuations l3ll. human motor acti vity |33 | and 
gait [33, HJ, IsE ISa, I33, meteorolo gy ISSl 1391 . c limate 
temperature fluctuations U El @ l43l l44 l4p . river 
flow and discharge [H |43 , electric signals [43, H^, ISOl , 
stellar x-ray binary systems 1^, neural receptors in bi- 
ological systems IsaKmusic [53 , and economics [5J, [S^, 
ISo. I57I IsM I59I l60l l6ll |. In many of these applications 
the main problem is to differentiate scaling features in 
a system's output which are inherent to the underlying 
dynamics, from the scaling features which are an artifact 
of nonstationarities or different types of transforms and 



filters. 

In two previous studies we have examined how different 
types of nonstationarities such as superposed sinusoidal 
and power-law trends, random spikes, cut-out segments, 
and patches with different local behavior affect the long- 
range correlation properties of signals [WL l62 | . Here we 
use the DFA method to investigate how the scaling prop- 
erties of noisy correlated signals change under linear and 
nonlinear transforms. Further, (i) we test to see under 
what types of transforms (filters) it is possible to derive 
information about the scaling properties of the signal of 
interest before the transformation, provided we know the 
correlation behavior of the transformed (filtered) signal, 
and (ii) we probe the "apparent" scaling of three common 
transformation functions after applying the DFA method 
— exponential, logarithmic and polynomial. We also 
evaluate the limitations of the DFA method under lin- 
ear and nonlinear transforms. Specifically, we consider 
the following: 

(1) Correlation properties of signals after transforms of 
the type: {xi} =^ {f{xi)}, where {xi} is a stationary 
signal with a priori known correlation properties. 

(i) Linear transform: {xi} =^ {axi + b}. Transforms 
of this type are often encountered in physical systems. 
For example: (a) from the fluctuations in the accelera- 
tion of a particle (measurable quantity), one can derive 
information about how the force (quantity of interest) 
acting on this particle changes in time without directly 
measuring the force: {a{ti)} => {F{ti) = ma{ti)}; 
(b) in pnp-transistors a difficult to directly measure 
base (input) current Ib (quantity of interest) is am- 
plified hundreds of times, so that small fluctuations in 
Ib may lead to significant (and measurable) changes in 
the collector (output) signal Ic (measurable quantity): 
{Ic{ti)} =^ {iBiU) = Ic{ti)/l3}, and (c) changes in the 
volume V (quantity of interest) of an ideal gas can be 
determined from fluctuations in the temperature (mea- 
surable quantity) provided the pressure is kept constant: 

{r(iO} =^ {vit.) = ^T(i,)}. 

(ii) Nonlinear polynomial transform: {xi} =^ {axf}, 
where k ^ 1 and takes on positive integer values. For 
example: (a) from fluctuations in the current / (mea- 
surable quantity) one can extract information about the 
behavior of the power lost as heat P (quantity of inter- 
est) in a resistor: {/(i,)} =^ {P{ti) = Rl^iU)}; (b) 
measuring the temperature T fluctuations of a radiating 
body the Stefan's law defines the power emitted per unit 
area: {Ti} =^ {ci = crT^}. Further, linear and nonlinear 
polynomial filters are also used to renormalize data series 
representing an identical quantity measured in different 
systems before performing correlation analysis, e.g., (i) 
normalizing heart rate recordings from different subjects 
to zero mean and unit standard deviation (linear filters), 
or (ii) extracting the absolute value (nonlinear filter) of 
the heartbeat fluctuations in datasets obtained from dif- 
ferent subjects p5| . 

In this study we consider two examples of nonlinear 
polynomial filters — quadratic and cubic filters — which 



represent the class of polynomial filters with even and 
odd powers, and we investigate how these filters change 
the correlation properties of signals. Since polynomial 
filters with even power wipe out the sign information in 
a signal, we expect quadratic and cubic filters to have 
a different effect. A recent study by Y. Ashkenazy et 
al. I25| shows that the magnitude of a signal (without 
sign information) exhibits different correlation properties 
from that of the original signal. Thus it is necessary to 
investigate how quadratic and cubic filters change the 
scaling properties of correlated signals. 

(iii) Logarithmic filter: {xi\ =^ {log(a;i + A)}, is also 
widely used in renormalizing datasets obtained from dif- 
ferent sources before comparative analysis. For exam- 
ple, to compare the dynamics of price fluctuations X{i) 
of different company stocks, which may have a different 
average price, one often first obtains the relative price 
returns R{i) — \oe.\X(i + 1)/X{i)] before performing cor- 
relation analysis |55l l63l | . It is assumed that upon taking 
the returns one does not alter the information contained 
in the original signal. To test this assumption we com- 
pare the correlation properties of the signal before and 
after a logarithmic filter. 

(2) Correlation properties of transformation functions 

When analyzing the correlation properties of a signal 
after a given transform, it may be valuable to know what 
is the DFA result for the transformation function itself. 
In addition, it is often the case that noisy signals are su- 
perposed on trends which can be approximated by a cer- 
tain function. Previous studies have demonstrated that 
the DFA result of a correlated signal with a superposed 
trend is a superposition of the DFA result for the signal 
and the DFA result for the analytic function represent- 
ing the trend [ly, IS^I ■ Here we investigate separately the 
results of the DFA for three functions which are very of- 
ten encountered in physical and biological processes: (i) 
exponential, (ii) logarithmic and (iii) power law. 

The layout of this paper is as follows: In Sec. m we 
describe how we generate signals with desired long-range 
power-law correlations and introduce the DFA method 
used to quantify correlations in nonstationary signals. In 
Sec. mil we compare the correlation and scaling proper- 
ties of signals before and after linear and nonlinear poly- 
nomial transforms. In Sec. IIVL we consider the effect of 
nonlinear logarithmic filter on the long-range correlation 
properties of stationary signals. In Sec. we investi- 
gate the performance of the DFA method on three an- 
alytic functions — exponential, logarithmic, and power- 
law — which are often encountered as trends in physical 
and biological time series. We systematically examine 
the crossovers in the scaling behavior of correlated sig- 
nals resulting from the transforms and trends discussed in 
Secs. lIIliVl the conditions of existence of these crossovers 
and their typical characteristics. We summarize our find- 
ings in Sec. I VII 



II. METHODS 

We analyze two types of signals: 

(1) stochastic stationary signals {2;^} {i = 
1, 2, 3, ...jiVmax) with different type of correlations 
(uncorrelated, correlated, and anti-correlated) and 
surrogate signals obtained from {xi} after linear and 
nonlinear transforms. We use an algorithm based on 
the Fourier transform to generate stationary signals 
jxj} with long-range power-law correlations as described 
in p2L |6J, |63- The generated signals {xi} have zero 
mean and unit standard deviation. 

(2) Exponential, logarithmic, and power-law functions 
which often represent transformations or trends in phys- 
ical and biological data. 

We use the detrended fluctuation analysis (DFA) 
method |a, l3| to quantify the correlation and scaling 
properties of these signals. The DFA method is described 
in detail elsewhere [13, |63| • Briefly, it involves the fol- 
lowing steps: (i) we integrate the signal after subtracting 
the global average; (ii) we then divide the time series 
into boxes of length n and perform, in each box, a least- 
square polynomial fit of order £ to the integrated signal 
to remove the local trend in each box; (iii) in each box 
we calculate the root-mean-square fluctuation function 
F(n) quantifying the fluctuations of the integrated sig- 
nal along the local trend; (iv) we repeat this procedure 
for different box sizes (time scales) n. 

A power-law relation between the average root-mean- 
square fluctuation function F(n) and the box size n in- 
dicates the presence of scaling: F{n) ~ n". The scale n 
for which this scaling holds represents the length of the 
correlation. The fluctuations in a signal can be charac- 
terized by the scaling exponent a, a self-similarity pa- 
rameter which quantifies the strength of the long-range 
power-law correlations in the signal. If a = 0.5, there 
is no correlation and the signal is uncorrelated (white 
noise); if a < 0.5, the signal is anti-correlated; if a > 0.5, 
the signal is correlated. Since we use a polynomial fit 
of order i, we denote the algorithm as DFA-£. Further, 
we note that for stationary signals {xi} with long-range 
power-law correlations, the value of the scaling exponent 
a is related to the exponent /3 in the power spectrum 
S{f) = f-P of signals {xj by /? = 2a - 1 0. Since the 
power spectrum is the Fourier transform of the autocor- 
relation function, one can find the following relationship 
between the autocorrelation exponent 7 and the power 
spectrum exponent (3: 7=1-/3 = 2 — 2a, where 7 is 
defined by the autocorrelation function C{t) — t~'^ and 
should satisfy < 7 < 1 . 

The upper threshold for the value of the scaling ex- 
ponent a is related to the order t of the DFA method: 
a < ^ -f- 1 for DFA-£ [T^. In addition, integrating the 
signal before applying the DFA method will increase the 
value of the scaling exponent a by 1, thus the upper 
threshold will become a-|-l < t+1 for DFA-£. Therefore, 
after integrating correlated signals with the scaling expo- 
nent a > t, one needs to apply the DFA method with an 



order of polynomial fit higher than £. We also note that 
for anti-correlated signals, the scaling exponent obtained 
from the DFA-€ method overestimates the true correla- 
tions at small scales [l^. To avoid this problem, one 
needs first to integrate the original anti-correlated sig- 
nal and then to apply the DFA-^ method [l3,E2- The 
correct scaling exponent a can then be obtained from 
F{n)/n [instead of F{n)] 10, 25, 62] . For that reason 
we first integrate and then apply the DFA method when 
considering anti-correlated signals. 



III. EFFECTS OF LINEAR AND NONLINEAR 
POLYNOMIAL TRANSFORMS 

In this section, we study the effect of linear and nonlin- 
ear polynomial transforms (filters) on the scaling proper- 
ties of stationary signals {xi} with long-range power-law 
correlations. Specifically, we consider two types of non- 
linear transforms — quadratic and cubic — as an exam- 
ple of even and odd polynomial filters. We generate the 
signals {xi} with linear fractal properties and with a pri- 
ori build-in correlations characterized by a DFA scaling 
exponent a |j, llOl |62| . We compare how the exponent a 
changes after the transform. 

We first test to see if these transforms affect the prop- 
erties of uncorrelated signals (white noise) . We find that 
the linear, quadratic, and cubic filters do not change the 
scaling properties of white noise — the curves of the de- 
trended fluctuation function F{n) for the flltered signals 
{f{xi)} collapse on the scaling curve of the original sig- 
nal {xi}, and the scaling exponent a = 0.5 remains un- 
changed [Fig.nia)]. 

For signals with correlations we find that the linear 
and nonlinear polynomial filters have a different effect. 
In particular, for both correlated (a > 0.5) and anti- 
correlated (a < 0.5) signals {xi} we find that the scaHng 
properties remain unchanged after the linear filter. In 
contrast, the quadratic and cubic filters change the scal- 
ing behavior of both correlated and anti-correlated sig- 
nals [Fig.^b-d)]. Specifically, for anti-correlated signals ^ 
we find that: (i) after the quadratic filter the scaling 
behavior is dramatically changed to uncorrelated (ran- 
dom) behavior with a = 0.5 at all scales; (ii) after the 
cubic filter the scaling (correlation) function F{n) of anti- 
correlated signals is also changed and exhibits a crossover 
from anti-correlated behavior at small scales to uncorre- 
lated behavior at larger scales [Fig.^b)]. We note that 
the quadratic filter removes the sign information in a sig- 
nal, thus completely eliminating the anti-correlations in 
a signal. In contrast, the effect of the cubic filter is not 
as strong as the effect of the quadratic filters, since a 
cubic filter preserves the sign information and the anti- 
correlations at small scales. For correlated signals we find 
that after both quadratic and cubic filters, the scaling be- 
havior is unchanged at small and intermediate scales. At 
large scales we observe a crossover to weaker correlations 
which is less pronounced when increasing the strength of 
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(a) White noise: a=0.5 

{2x^+10} 

{x,^} 

-— {Xi=} 




10' 



10" 



a=0.5 



s 



DFA-1 



10" 



10' 



10' 



10^ 



10' 



10" 



10"' 



(b) Anti-correlated signals: a=0.2 
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FIG. 1; Effects of linear, quadratic, and cubic filters on the scaling behavior of long-range correlated stationary, Gaussian 
distributed (zero mean and unit standard deviation) signals {xi}: (a) uncorrelated, (b) anti-correlated, (c) correlated, and (d) 
strongly correlated. The length of each signal is Nmax = 2^^. In our analysis we use the DFA-1 method. The curves of the 
detrended fluctuation function F{n) for all signals are vertically shifted so that they start at the same value of F{n) at the 
minimal scale n. For anti-correlated signals we first integrate and then apply the DFA-1 method to avoid overestimation of the 
true correlations at small scales due to limitations of the DFA method (|il0l.62l and Sec. [HJ. Our analysis shows that after a 
linear filter the scaling behavior remains unchanged. In contrast, nonlinear polynomial filters change the scaling behavior of 
anti-correlated and correlated signals, leading to crossovers at large scales. 



the correlations (higher values of a) in the signal {xi} 
[Fig. ^c-d)] . For signals with very strong correlations 
(a > 1), we find that the scaling behavior remains al- 
most unchanged after nonlinear polynomial filters. We 
also find that the quadratic filter leads to a more pro- 
nounced crossover at large scales compared to the cubic 
filter for all positively correlated signals. 



IV. LOGARITHMIC FILTER 

In addition to nonlinear polynomial transforms, loga- 
rithmic transforms are often used in preprocessing proce- 
dures when there is a need to renormalize output signals 
obtained from different systems before comparing their 
correlation properties |55J. In this section, we investigate 



the effect of logarithmic filters on the scaling properties 
of stationary signals with long-range power-law correla- 
tions. 

We first generate stationary correlated signals {xi} 
with a zero mean and unit standard deviation, and with a 
priori known and controlled correlation properties quan- 
tified by DFA scaling exponent a. To ensure that all 
values in the signal are positive, before the logarith- 
mic transform, we shift {xi} =^ {xi + A}, where A = 
—Xmin + e, Xrain IS the minimal value in the series {xi\ 
and e is a positive constant. This linear transform does 
not alter the correlation properties of {xi}, as demon- 
strated in Sec. IIIII Fig. ^ Next we integrate the signal 
after the logarithmic transform {log^gCa^i— a^min + e)} and 
we perform a DFA-2 analysis. 

For uncorrelated (white noise) signals after the loga- 
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FIG. 2: Effects of the nonlinear logarithmic filter {log^Q{xi — Xmin + e)} on the scaling behavior of stationary correlated signals 
{xi}, where Xmin is the minimal value in the original signal {xi} and e is a positive constant. The original signals {xi} have zero 
mean, unit standard deviation, and length Nmax ~ 2^**. (a) Original strongly anti-correlated signal {xi} with DFA correlation 
exponent a = 0.1 and the corresponding signal after logarithmic filter, (b) DFA scaling curves F{n) for anti-correlated signals 
and white noise after the logarithmic filter show a crossover to "white noise" behavior (i.e., slope= 0.5) at large scales. To 
obtain more accurate scaling, we first integrate the signal {\og-^Q{xi — Xmin + ()} and then apply DFA-2 method (see Sec.^. (c) 
Original strongly correlated signal {xi\ with the DFA correlation exponent a = 2 and the corresponding signal after logarithmic 
filter, (d) DFA scaling curves for correlated signals {xi} after the logarithmic filter. We find that the logarithmic filter does 
not change the correlation properties of signals with certain positive correlations (e.g., a — 0.7 and a = 1.2), though it weakens 
the correlations in signals with very strong positive correlations (e.g., a = 2). 
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(b) Correlated signals: a=2 
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FIG. 3: Dependence of the effect of logarithmic filter {\ogii~,{xi — Xmin + e)} on the offset parameter e. (a) Detrended 
fiuctuation function F{n) from the DFA-2 after integration of {logjg(a::i — Xmin + e)}, for an anti-correlated signal with the 
DFA correlation exponent a = 0.1 and varied values of e. We find that for smaller values of e, there is a more pronounced 
crossover to uncorrelated behavior with a = 0.5. (b) Detrended fluctuation function F{n) from the DFA-2 after integration of 
{\ogiQ{xi — Xmin + t)} , for a signal with strong positive correlations (a = 2) and varied values of e. We find that signals with 
strong positive correlations are less affected by the logarithmic filter compared to the anti-correlated signals in (a) and that for 
smaller values of e, there is a more pronounced crossover. 
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FIG. 4: Relation between the scaling exponent a of the orig- 
inal "input" stationary signals and the correlation exponent 
aout of the signals after the logarithmic filter {log(xi — Xmin + 
e)}, where Xmin is the minimal value in the original signal {xi} 
and e is a positive constant. Uout is obtained from the DFA-2 
analysis after integrating the signal {log(xi — Xmin + e)} and 
fitting the detrended fluctuation function F{n) in the region 
n £ [30, 3000]. Our results show that for signals with a corre- 
lation exponent a outside the shaded region, the logarithmic 
filter changes the scaling behavior (aout i= a) and this change 
depends on the offset parameter e. 



rithmic filter, we find no change in the scaling proper- 
ties and the correlation exponent remains a = 0.5 in 
the entire range of scales [Fig. [2Ib)]. However, we find 
that the scaling properties of signals with certain de- 
gree of correlation change significantly. Specifically, for 
anti-correlated signals (a < 0.5) we observe a crossover 
to uncorrelated (white noise) behavior at large scales. 
This crossover becomes more pronounced (and shifted 
to smaller scales) when increasing the strength of anti- 
correlations (decreasing a) [Fig. |2Ib)]. This crossover 
behavior is caused by negative spikes in the signal fol- 
lowing the logarithmic transform [Fig. |2Ia)]. A similar 
effect was previously reported for stationary correlated 
signals with superposed random spikes |62|- For corre- 
lated signals (a > 0.5), we find a threshold value for the 
correlation exponent (Xth ~ 1-3, below which the scal- 
ing properties of the signal remain unchanged after the 
logarithmic filter. Above ath there is a reduction in the 
strength of the positive correlations, i.e., the value of the 
estimated exponent after the logarithmic filter is much 
lower compared to the correlation exponent a in the orig- 
inal signal [Fig. |2Id)] . 

Since the logarithmic filter is a nonlinear transform 
which diverges for values of the signal {xi — Xmin + e} 
close to zero, we next test how the scaling properties 
of the signal depend on the value of the offset parame- 
ter e. We consider anti-correlated and correlated signals 
with fixed values of a and varied e. For strongly anti- 
correlated signals we find that even for large values of e. 



there is a crossover to uncorrelated behavior in the scal- 
ing curve F{n) at large scales (note that e is the minimal 
value of the signal {xi — Xmin + e})- This crossover shifts 
to smaller scales with decreasing e [Fig|3fa)]. Further, we 
find that for decreasing e, the scaling curves F{n) con- 
verge to a single curve, indicating random uncorrelated 
behavior in the range of large and intermediate scales. 
For anti-correlated signals with a = 0.1 we find that this 
convergence is reached for e < 0.1 [FiglSJa)]. For sig- 
nals with strong positive correlations [a > ath), we also 
observe a change in the scaling behavior which becomes 
more pronounced when e decreases. However, in contrast 
to the anti-correlated signals, the deviation from the ex- 
pected accurate scaling starts at intermediate scales and 
extends to smaller scales with decreasing e [Fig|3fb)]. For 
signals with very strong correlations, e.g. a = 2, the 
deviation from the accurate scaling is observed only for 
e < 0.1, while for e > 0.1, there is no effect on the scaling 
[FigEIb)]. This is in contrast to the situation observed 
for signals with strong anti-correlations (a = 0.1) where 
the logarithmic filter alters the scaling behavior even for 
much larger values e > 10 [FiglSJ^a)]. 

Finally we study the relation between the scaling ex- 
ponent a of the original "input" signal and the estimated 
exponent aout of the "output" signal after the logarith- 
mic filter. We find that for correlated signals within given 
range for the value of the scaling exponent a G [0.4, 1.3), 
there is no change in the scaling properties after the log- 
arithmic transform. However, for signals with correla- 
tion exponents a < 0.4 and a > 1.3, we find that the 
logarithmic transform can dramatically change the scal- 
ing behavior and this effect also strongly depends on the 
value of the offset parameter e [Fig0]. Therefore, the 
logarithmic filter is not recommended for anti-correlated 
signals and signals with very strong positive correlations 
— applying this filter will mask the true correlations in 
the original signals. 



V. RESULTS OF THE DFA FOR 
TRANSFORMATION FUNCTIONS 

In this section we investigate the scaling properties of 
three functions: exponential, logarithmic, and power-law. 
These functions are often used in signal processing as 
transforms of various stochastic correlated signals and 
also appear as trends superposed on noisy signals de- 
rived from physical and biological systems. In previous 
work 1101.1631 we have demonstrated that the scaling be- 
havior of a correlated signal with a superposed trend is 
superposition of the scaling behavior of the correlated sig- 
nal and the "apparent" scaling behavior obtained from 
the DFA method for the analytic function representing 
the trend. Therefore, understanding the results of the 
DFA for certain analytic functions becomes a necessary 
step to quantify the scaling behavior of system's outputs 
where correlated fiuctuations are superposed with differ- 
ent trends. 
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FIG. 5: The results of the DFA method for general exponential function: y = exp(ca; + a), < a; < l,a; = i/Nmax,i = 
1, 2, ..., Nmax,Nm,ax = 2^^, where c = ±1 and offset a is a positive constant, (a) Detrended fluctuation function F{n) obtained 
using the DFA-2 method for different values of the offset parameter a. While there is a vertical shift in F{n) for different values 
of a, all scaling curves are characterized by an identical slope a. (b) Dependence of the scaling exponent a on the parameters 
a and c. We find that for any exponential function the scaling exponent ct depends only on the order (. of the DFA method: 
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FIG. 6: The results of the DFA method for general logarithmic function y — logj^Q(x + a),0 < x < l,x — i/Nmax,i ~ 
1,2, ...,Nmax,Nmax = 2^"^, whcrc offset a is a positive constant, (a) Detrended fiuctuation function F{n) obtained using the 
DFA-2 method for different values of the offset parameter a. We find that the slope of the scaling curve (scaling exponent 
a) depends on the value of the offset a. (b) Dependence of the scaling exponent a on the offset a [fitting region for a is 
n £ (30, 3000)]. We observe a dramatic change from a = 1.5 ato~0toa = i-|-lata> 0.01, where £ is the order of the DFA 
method. 



(i) We first consider the exponential function in the 
form: y — exp(ca; -|- a), where < a; < 1, a; = i/Nmax, 
i = 1, ...,Nraax, Nmax = 2^''', the parameter c = ±1, the 
offset parameter a is a positive constant. We show the 
result of the DFA method in Fig. |5| We find that the 
slope of the detrended fluctuation function F{n) vs. the 
scale n obtained from the DFA method does not depend 
on the values of the parameters c and a (there is only 
a vertical shift in F{n) for different values of a and c) 
[Fig.jSJa)]. Instead, we find that the DFA scaling expo- 



nent a depends only on the order £ of polynomial fit in 
the DFA method — a = £+1 — suggesting that the re- 
sults of the DFA method do not depend on the details of 
the exponential function [Fig.[SJb)]. An analytic deriva- 
tion for the fluctuation function F(n) and the value of the 
scaling exponent a obtained from DFA-1 is presented in 
Appendix A. 

(ii) We next consider the performance of the DFA 
method on a logarithmic function of the general form: 
y = logio(a; + a), where < a; < 1, a; = i/N^ax, 
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FIG. 7: The results of the DFA method for general power-law function y — {x + a) ,Q < x < l,x = i/Nmax,i ~ 
1,2, ..., Nmax, Nmax = 2^^, where A is the power and the offset parameter a is a positive constant, (a) Detrended fluctua- 
tion function F{n) obtained using the DFA-2 method for fixed A — —0.39 and different values of the offset parameter a. We 
find that the slope of the scaling curve (scaling exponent a) depends on the value of a. (b) Dependence of the scaling exponent 
a on the offset a for different values of the power A [fitting region for a is n £ (30, 3000)]. We observe a dramatic change from 
Of ~ A + f .5 at a ~ for different values of A to a = ^ + 1 for a > fO~^, where £ is the order of the DFA method. 



i — 1, ..., Nfnax, Nmax = 2"'^^ and the offset parameter 
a is a positive constant. Specifically, we investigate the 
dependence of the DFA scaling exponent a on the value 
of the offset parameter a. We find that for very small 
values of a, the DFA scaling exponent is a — 1.5. With 
increasing a, we observe a crossover in F{n) at interme- 
diate scales n — from a = 1.5 at large scales to a = 3 at 
small scales for DFA-2 [Fig.|SJa)]. For larger values of a, 
we observe a scaling behavior in F(n) characterized by 
a single exponent a = 3 in the entire range of scales n 
[Fig- EI a-)] • In Fig-EIb) we present the dependence of the 
DFA scaling exponent a [obtained in the fitting range 
n G (30, 3000)] on the offset parameter a for different 
DFA order £. We find that for a < 10~^ the exponent 
a does not depend on the order £ of the DFA method 
and takes on a single value a = 1.5. In contrast, for 
large values of a > 10~^, the exponent a depends only 
on the order £ of the DFA method and takes on values 
a = £ + 1. This behavior is identical with the behavior 
obtained for the exponential function in Fig. [SJb). For 
intermediate values of a, we observe a crossover in the 
scaling behavior of the fluctuation function F{n) from 
a = 1.5 to a = ^+ 1. 

(iii) Finally, we consider the general power-law func- 



tion: y = (x -\- a) , where < a; < 1, x = i/N^ax, 
i — I, ■■■,Nmax, Nmax = 2^^, the powcr A takes on real 
values and the offset parameter a is a positive constant. 
As in the case of the logarithmic function, we find again 
that the DFA scaling exponent a depends on the value 
of the offset parameter a [Fig. Oa)] . For certain fixed 
values of A and with increasing a, we observe a grad- 
ual transition in the fiuctuation function F{n) from a 
scaling behavior spanning over a broad range of scales 
n characterized by a small value of the exponent a to a 



crossover at intermediate scales n for larger values of a, 
and finally to a scaling spanning over all scales n with 
exponent a = 3 for large values of a for DFA-2. In a 
previous study [l^] we have found a specific relationship 
between the DFA exponent a and the value of the power 
A for the case of power-law function with offset param- 
eter a = 0: a^£+liorX>£- 0.5; a ~ A -I- 1.5 
for -1.5 < A < ^ - 0.5; a = for A < -1.5, where 
£ is the order of polynomial fit in the DFA-^ method. 
Our current analysis shows that this behavior is even 
more complicated when a > [Fig. [7|[b)]. Specifically, 
we find that for values of a < 10^^ the scaling exponent 
a [obtained in the fitting range n e (30, 3000)] depends 
only on the value of the power A: a ~ A -I- 1.5. In con- 
trast, for large values of the offset parameter a > 10~^, 
we find that the exponent a depends only on the order 
£ of the DFA method, and takes on values a — £ + 1, 
which is similar to the results obtained for the general 
exponential and logarithmic functions in this range of a 
[Fig. Efb) and Fig. Elb)]). For intermediate values of a 
and for —1.5 < A < £ — 0.5, we observe a crossover in the 
scaling behavior of the fiuctuation function F{n) from 
a ~ X -\- 1.5 to a — £ -h 1- Further, we find that for 
X > £ — 0.5, the DFA-^ scaling exponent remains con- 
stant a = £ + 1 , and does not depend on the values of the 
offset parameter a — we note that for A = 0.41 (close to 
A = 0.5 = £ — 0.5 for DFA-1) the dependence of a on a 
is close to a horizontal line [Fig.[7Ib)]. 

Analytic arguments 

Our results show that for large values of the offset pa- 
rameter a, the detrended fluctuation function F{n) for 
all three analytic functions — exponential, logarithmic, 
and power-law — exhibits identical slope, where the DFA 
scaling exponent a does not depend on the particular 



functional form but only the order i of the DFA method: 
a = ^ + 1 [Fig. inib), Elb) andHb)]. The reason for 
this common behavior is that (i) for large values of a, 
in each DFA box of a given length n, all three functions 
can be expanded in converging Taylor series, allowing for 
a perfect fit by a finite order polynomial function, and 
(ii) that, due to this convergence, the same polynomial 
function can be used when shrinking the box length n. 
In contrast, for very small values of the offset parameter 
a, the DFA results for all three functions are distinctly 
different and does not depend on the order I of the DFA 
method. Below we give some general analytic arguments 
for the dependence of the DFA exponent a on the offset 
parameter a presented in Figs. El and 

(i) General exponential function: y — exp(x + a),0 < x < 
1. 

First, we substitute the variable x hy z = x + a: y = 
e^, z e (a, 1 + a]. Next, we consider a DFA box starting 
at the coordinate z' = s and ending at z" = s + 1, where 
t is proportional to the number of points n in the box 
— i = (1 + a - a)n/Nmax = n/N,nax- For any value 
of z G (s, s + t) we can expand the function in a Taylor 
series: 



exp(s + zo)|o<^o<t = ^' 



"°+2! 



(1) 



Since this expansion converges, a finite polynomial 
function can accurately approximate the exponential 
function in each DFA box. We note that the DFA-^ 
method applied to above polynomial functions gives the 
scaling exponent a = £ + I (see [l3)- Thus, for any 
exponential function we find that the DFA scaling does 
not depend on the value of the offset parameter a and 
depends only on the order £ of the polynomial fit in the 
DFA-£ procedure [Fig.E^b)]. 

(ii) General logarithmic function: y — log]^o(^ + a),0 < 
X <1. 

First, we substitute the variable x hy z — x + a: y ~ 
logiQ{z),z e (a, 1 + a]. Next, we consider a DFA box 
starting at the coordinate z' — s and ending at z" — s+t, 
where t is proportional to the number of points n in the 
box — t = n/Nmax- For any value of z £ {s,s + 1) the 
Taylor expansion is: 



login (z) = 



logio(s + ^o)|o<zo<t 
ln(l + zo/s) 

s 2 



& 



(?)' 



•(2) 



This series is converging only when zq/s < 1, i.e., zq < 
s. Since zq G (0, t), the condition for convergence in any 
DFA box (s, s + t) partitioning the function is i < s. 
From t = n/Nmax and s G [a, 1 + a — t], we find that if 
a > n/Nmax, the logarithmic function in all DFA boxes 
is converging, and thus each box can be approximated 



by a polynomial function, leading to scaling exponent 
a — £-\-\ — depending only on the order I of the DFA-^ 
method [Fig.|S]. 

When t > s, for certain values of zq G (0, f), the series 
in Eq. Q is diverging. Since s G [a, 1 + a — t], for s = 
a < t = n/Nmax, we find that the logarithmic function 
is divergent in the first DFA box {a, a + t), leading to 
deviation in the DFA scaling for small values of a [Fig.|S]. 

(iii) General power-law function: y = {x + a)'^, x G (0, 1]. 
First, we substitute the variable a;byz = a; + a: y — 
z'^, z G (a, 1 + a]. Next, we consider a DFA box starting 
at the coordinate z' = s and ending at z" = s + 1, where 
t is proportional to the number of points n in the box — 
t = n/Nmax- For any value of z G (s, s + t) the Taylor 
expansion is: 



zo^ ^ 
s 



1 



= 1 + A 



£o A(A - 1 ) f zq 
s 2! 



(3) 



Similar to the case of the logarithmic function, this 
series is converging only when zq/s < 1. Following the 
same arguments as for the logarithmic function we find 
that when a > n/Nmax, the power-law function is con- 
verging in any DFA box, and thus can be approximated 
by a polynomial function, leading to the scaling expo- 
nent a ^ £ + I [Fig.[7|, which is identical to the case of 
exponential and logarithmic function. 

In contrast, for a < n/Nmax, the power-law function 
is divergent in the first DFA box (a, a + t), as in the case 
of the logarithmic function, leading to a deviation in the 
scaling of F{n) for small values of a [Fig. [7]. While in 
the case of logarithmic function this divergence leads to 
a fixed scaling exponent a = 1.5, for power-law functions 
the value of the scaling exponent a depends also on the 
power A [Fig.!?]. 

We note that the above arguments can be used to esti- 
mate the results of the DFA method for other functions. 
For all functions which can be expanded in convergent 
Taylor expansion of a polynomial form in each DFA box 
partitioning the function, the DFA method leads to iden- 
tical scaling results with the exponent a = ^ -I- 1 , which 
is a notable inherent limitation of the method. When 
there is divergent behavior in some or all of the DFA 
boxes partitioning a function, the DFA scaling exhibits 
crossover behavior to different values of the scaling ex- 
ponent a which depends on the functional form and the 
specific parameters of the function. 



VI. CONCLUSIONS 

In summary, our study shows that linear transforms do 
not change the scaling properties of a signal. However, 
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the correlation properties of a signal change after apply- 
ing a polynomial filter. Moreover, such change depends 
on the type of correlations (positive or anti-correlations) 
in the signal, as well as on the power (odd or even) of the 
polynomial filter. For the logarithmic filter we find that 
the scaling behavior of the transformed signal remains 
unchanged only when the original signal satisfies certain 
type of correlations (characterized by scaling exponent 
within a given range). Comparing the "apparent" scal- 
ing behavior of the exponential, logarithmic, and power- 
law functions we find that within certain range for the 
values of the parameters, the DFA fluctuation function 
F{n) exhibits an identical slope, and that the DFA re- 
sults of a class of other analytic functions can be reduced 
to these three cases. We attribute this behavior to spe- 
cific limitations of the DFA method. Therefore, careful 
tests are necessary to accurately estimate the correlation 
properties of signals after nonlinear transforms. 



values. The first step of the DFA method is to integrate 
the signal [SecHT]: 



exp(— +a)c^2/ = iV- 



(Al) 



where N is the length of the signal and x G (0,iV]. 
We divide the variable in the exponential by A^, so that 
(x/N) is in the interval (0, 1], as considered in Sec. 
The next step of the DFA method is to divide the in- 
tegrated signal into boxes of length n. For DFA-1, the 
squared detrended fluctuation function in the fc— th box, 
F^{n,k), is 
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APPENDIX A: DFA-1 IN EXPONENTIAL 
FUNCTIONS 

We consider an exponential function of the type 
exp{cx + a), where the parameters c and a take on real 



1 



F^{n,k) = - 

" ■l(k-l)n 



kn 



N- 



[bkX - dk) 



dx, 
(A2) 



where the parameters b^ and dk are obtained by a lin- 
ear fit to the integrated signal using least squares in the 
fc— th box. These two parameters can be obtained analyt- 
ically, although their expressions are too long. To obtain 
the squared detrended fluctuation function for the entire 
signal partitioned in non-overlapping boxes of length n, 
we sum over all boxes and calculate the average value: 



J 



T N/n N/n .„ 



N- 



[bkX - dk) 



dx. 



(A3) 



Here, the index k in the sum ranges from 1 to N/n 
(there are N/n boxes of length n in the signal of length 
N) . Using the analytical expressions for hk and dk , F^ (n) 
can be presented analytically in the form: 



where 



F'^{n) ^g{n) ■ h{n), 



(A4) 



lin) = {- 



g{n) = \ -SNc^n'^ (l + e"^ + e^) + c^n^ (e^ ~ ^) + ^"^^^ 



-(e"-l) N - cn + cne " 



■]} (A5) 



r 



and 



Due to the complexity oi g(n) and h(n), the expression 



h{n) 



g2a (g2c _ 1) ^2 

2c6 f e^r - 1 j n3 



(A6) 



of F'^{n) is very complicated. However, as n < TV (and 
usually, n <C iV), one can expand F'^{n) in powers of n 
to obtain: 



F\n) 



c (e2- - 1) e'"' 

— -^ 7 

14407V2 



(A7) 



Finally, for the detrended fluctuation function F{n) we 
obtain: 
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Thus the DFA-1 scaling exponent is a = 2 (in agree- 
ment with the numerical simulation in Sec. [3 Fig.|SJl. In 
general, we can obtain in a similar way that a — £ + I, 
when DFA-^ with an order £ of polynomial fit is used. 



^t X /c(e2=-l)e° 2 

F(n) ~ \ — — n^. 

^ ' ^' 1440 N 



(A8) 
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